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By an extension of the Feynman-Kleinert variational approach, we calculate the temperature- 
dependent effective classical potential governing the quantum statistical properties of a hydrogen 
atom in a uniform magnetic field. In the zero-temperature limit, we obtain ground state energies 
which are accurate for all magnetic field strengths from weak to strong fields. 



I. INTRODUCTION 

o . 
o- 

The recent discovery of magnetars has renewed interest in the behavior of charged particle systems in the presence 
CS| \ of extremely strong external magnetic fields jjj. In this new type of neutron stars, electrons and protons from decaying 
neutrons produce magnetic fields B reaching up to 10 15 G, much larger than those in neutron stars and white dwarfs, 
03 \ where B is of order 10 10 - 10 12 G and 10 6 - 10 8 G, respectively. 

Analytic treatments of the strong-field properties of an atomic system are difficult, even in the zero-temperature 
limit. The reason is the logarithmic asymptotic behavior of the ground state energy In the weak- field limit, 

perturbative approaches yield well-known series expansions in powers of B 2 . These are useful, however, only for 
B < B , where B is the atomic magnetic field strength B = e 3 M 2 /h 3 w 2.35 x 10 5 T = 2.35 x 10 9 G. 

So far, the most reliable values for strong uniform fields were obtained by numerical calculations An analytic 
mapping procedure was introduced in Ref. J3j to interpolate between the weak- and strong-field behavior, and a 
I variational approximation was given in Ref. |7[], both with quite good results. 

In this note, we use an extension of the Feynman-Kleinert variational approach |g[ to find a single analytic approx- 
t-H ' imation to the effective classical potential of the system for all temperatures and magnetic field strengths. From this, 
the quantum statistical partition function can be obtained by a simple configuration space integral over a classical- 
looking Boltzmann-factor. In the zero-temperature limit, the effective classical potential is the ground state energy 
' of the system. 

^: 

Ph. II. EFFECTIVE CLASSICAL POTENTIAL 

H— > ' 

The Hamiltonian of the electron in a hydrogen atom in the presence of a uniform external magnetic field pointing 
^ ' along the positive z-axis is 

>: ff (p,x) = ^-i^( P) x) + i^-^. (2.1) 

Here we have used the symmetric gauge A(x) = (B/2)(—y, x, 0), and denoted the z-component of the orbital angular 



momentum by l z (p,x) — (x x p) 2 . The quantum statistical partition function can always be expressed as a classical- 
looking configuration space integral f| 

d 3 x 



Z = I exp[-/3V eff (x )], (2.2) 



where Ath = y 27r?i 2 /3/Af is the thermal wavelength, j3 = 1/fcsT is the inverse temperature, and T4ff(xo) is the 

effective classical potential T4a(xo). Generalizing the development in Ref. B, this is defined by the phase space path 
integral 
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exp[-/3y cff (xo)] = X 3 th J d 3 Po Iv 3 xV 3 pS( Xo -x(r))d(p - p(r)) e~ A ^/ n , (2.3) 
where ^4[p,x] is the Euclidean action 

rhP 

A[vM= drHp(r)x(T)+fl(p(r),x(T))], (2.4) 
Jo 



and x(t) = g?tx(t)/?i/3 and p(r) = J Q drp(r)/?i/3 are the temporal averages of position and momentum. The 
special treatment of xo and po is necessary, since the classical harmonic fluctuation widths (x 2 ) cl and (p 2 ) cl are 
proportional to the temperature T (Dulong-Petit law). Thus they diverge for T — > oo and their fluctuations cannot 
be treated pertubatively. In contrast, the fluctuation widths ((x — xo) 2 ), ((p — Po) 2 ) around xo and po go to zero for 
large T and are limited down to T = 0, thus allowing for a treatment by variational perturbation theory For this 
we rewrite the action (2.4) as 

A[p, x] = A%°'*° [p, x] + Ant[ P , x], (2.5) 

with a harmonic trial action 

rh/3 



/|PO,Xo 



f np r l l 

P,x] = J dr|-i[p(r) - p ] • x(r) + — [p(r) - p ] 2 + -fi_Li/ z (p(r) - p , x(t) - x ) 

+ iMOi 2 [x» - X( [] 2 + iikTOf [*(r) - ^o] 2 }, (2.6) 

in which x- 1 = {x, y) denotes the transverse part of x. The frequencies fl = (Oj_x, f2j_2 5 are arbitrary for the mo- 
ment. Inserting the decomposition (2.5) into ( |2.3| ), we expand the exponential of the interaction, exp {— ^li nt [p, x]/?i}, 

(2.7) 



yielding a series of expectation values of powers of the interaction 

M£ lt [p,x])P- X0 = fv 3 xV 3 pA? nt [p,x}6(x -^r))6(p -p^))ex P |-i^' X0 [p,: 

The path integral over the Boltzmann-factor involving the harmonic action (2.6) is exactly solvable and yields the 
restricted partition function 

po , xo = hpn + /2 my2 wyi 

n smhhf3n+/2 sinhn/30_/2 sinh/^Qu/2' { '' 

where Q± = \Q±i ± f2j_2|/2. Rewriting the perturbation series as a cumulant expansion, evaluating the expectation 
values, and integrating out the momenta on the right-hand side of Eq. ( |2.3|) leads to a series representation for 
the effective classical potential VJ.fr (xo). Since it is impossible to sum up the series, the perturbation expansion 

must be truncated, leading to an .ZVth-order approximation Wfj (xo) for the effective classical potential. Since the 

parameters Q are arbitrary, W^'(xq) should depend minimally on Q. This determines the optimal values = 

(f2^(xo), (xo), fin (xo)) of Nth order. Reinserting these into Wf^\xo) yields the optimal approximation 

WW(x ) = lV s g ) (xo)._ 

The first-order approximation to the effective classical potential is 

V7( 1) (x )=-ilnZP- Xo + (c, c -f7 ±1 )fei(x )-i(f7i 2 - W 2 ) a i(x )-iAm 2 a 2 (x )-(^y°' ^ (2.9) 



where the quantities a^_(xo), ^(xo), and ajj(xo) are the transverse and longitudinal fluctuation widths 

al(^) = {x"(r)) P n °' x \ a 2 (x ) = (z 2 (T))^ Xo , 6 2 x (x ) = < x(r) Py (r) . (2.10) 



The expectation value of the Coulomb potential on the right-hand side of Eq. (2.9) has the integral representation 



7T II J a|(x ) + £ 2 [a 2 L (x ) - a 2 (x )] 

£2 / J i „,2 



N " M) ^ ^2 ( a 2 (x ) + £>i(xo) - a 2 (x )] + a 2 (x ) 1 ^ ' ^ ' ' ' 
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The variational energy (2.9) is minimized at each xo, and the resulting W( n \xq) is displayed for a low temperature 
and different magnetic fields in Fig. [I]. 

From now on we set Ti = e 2 = ks = c = M = 1. Thus, energies are measured in units of eo = Me 4 /7i 2 = 2Ryd ~ 
27.21 eV, temperatures in eo/fcs ~ 3.16 x 10 5 K, distances in Bohr radii olb = h /Me 2 s» 0.53 x 10~ 8 cm, and magnetic 
field strengths in B = e 3 M 2 /h 3 ps 2.35 x 10 5 T = 2.35 x 10 9 G. 
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FIG. 1. Effective classical potential plotted along along two directions: once as a function of the coordinate po = y/ Xq + Vq 
perpendicular to the field lines at zo = (solid curves), and once parallel to the magnetic field as a function of zo at po = 
(dashed curves). The inverse temperature is fixed at j3 = 100, and the strengths of the magnetic field B are varied (all in 
natural units). 



The curves W^'(xo) are plotted to show their anisotropy with respect to the magnetic field direction. The 
anisotropy grows when lowering the temperature and increasing the field strength. Far away from the proton at 
the origin, the potential becomes isotropic, due to the decreasing influence of the Coulomb interaction. Analyt ically , 
this is seen by going to the limits po — > oo or zq — > oo, where the expectation value of the Coulomb potential (2.11) 
tends to zero, leaving an effective classical potential 



--lnZg"*° + (u c 



2 II II 



(2.12) 



This is xo-independent, and optimization yields the constants = = u c and = 0, with the asymptotic 
energy 



^ (1) (xo) 



-4m 



hf3u c /2 



[3 sinh h/3uj c /2 ' 



(2.13) 



The B = -curves are, of course, identical with those obtained from variational perturbation theory for the hydrogen 
atom 0. 

For large temperatures, the anisotropy decreases since the violent thermal fluctuations have a smaller preference of 
the z-dircction. 



III. ZERO-TEMPERATURE LIMIT 



At zero-temperature, the first-order effective classical potential (2.9) at the origin yields an approximation for the 



ground state energy of the hydrogen atom in a uniform mag nctic field: Eg* = limp^ w£\o) 
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40 



"/ ~ A \ 1^1 



(3.1) 



with the expectation value for the Coulomb potential 



1 



■ arctan • 



±2 



X < 



"lb 



[ 2i Y n,| - n ±2 i - 2 ^20||/o ±2 - 1' 



20ii > r2j_2j 
2f2|| = S2j_2, 

2fi|| < fl ±2 - 



(3.2) 



Equations (3d) and (|3.2| ) are independent of the frequency parameter f2j_i, such that optimization of the ground state 



energy (3.1) is ensured by minimization. Reinserting the extremal f2j^ and Ql, 1 ' into Eq. (|3.l[) yields the first-order 



approximation to the ground state energy E^\B). In the absence of the Coulomb interaction the optimization with 



respect to f2j_ 2 yields Q 



(i) 

±2 



0J C 



rendering the ground state energy E^\B) = oj c /2, which is the zeroth Landau level 



in this special case. The trial frequency fin must be set equal to zero to preserve translational symmetry along the 
z-axis. 

In the opposite limit of a vanishing magnetic field, Eq. (3.1) coincides with the first-order variational result for 
the ground state eneray of the h ydr ogen atom, whose optimization gave E^\B = 0) = — 4/37T w —0.4244 [2 Ryd] 
obtained in Refs. MM. In Ref. |0L the B = -system was treated up to third order leading to the much more 
accurate result E^'(B = 0) « -0.490 [2Ryd], very near the exact value E ex (B = 0) = -0.5 [2Ryd]. 

Let us investigate the asymptotics in the strong- field limit B — ► oo. The ^-dependence of the binding energy 



e{B) = 



B 



E 



(3.3) 



is plotted in Fig. g, where it is compared with the results of Ref. ||, with satisfactory agreement. Our results 
are of similar accuracy as those of other first-order calculations, for example those from the operator optimization 
method in first order of Ref. . The advantage of variational perturbation theory is that it yields good results for all 
magnetic field strengths. From our experience with the fast convergence of the method ||, Chaps. 5,9], higher orders 
of variational perturbation theory will push the approximations rapidly towards the exact value. 
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FIG. 2. First-order variational result for the binding energy (3.3) as a function of the strength of the magnetic field. The 
dots indicate the values of Ref. The dashed curve shows the simple estimate 0.5 In 2 B of Landau-Lifschitz Q]. 
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A. Weak-Field Behavior 



The calculations of the binding energy for weak magnetic fields show that the ratio i] = 2f2||/f2j_2 is always smaller 
than one if B ^ 0. Setting O = f2j_2, we rewrite the binding energy as a function of Q and 77: 

:Wfm«! "fix? 




This is minimized with respect to the new variational parameters 77 and Q by expanding rj(B) and f2(B) in powers of 
B 2 with unknown coefficients, and inserting these expansions into extremality e qua tions. The expansion coefficients 
are then determined order by order. The optimal expansions are inserted into (3.4), yielding the optimized binding 
energy e^(B) as a power series 

00 

c {1) (B) = j -J2e n B^. (3.5) 

The coefficients e n are listed in Table Q and compared with the exact ones of Ref. [|. Of course, the higher-order 
coefficients of this first-order variational approximation become rapidly inaccurate, but the results can be improved, 
if desired, by going to higher orders in variational perturbation theory as in Ref. |9j, Chaps. 5,9]. 



B. Strong-Field Behavior 



In the discussion of the pure magnetic field we have mentioned that the variational calculation for the ground state 
energy, which is associated with the zeroth Landau level, yields a frequency fij^ oc B, while fiii = 0. We therefore 
use the assu mpti ons fij_ = il±2 ^S> 2f2ii and On <C B for an analytic study of t he strong-field behavior of the ground 



state energy (3.1). We expand the last expression of the expectation value fl3.2|) in terms of 2f2||/f2j_, and reinsert this 



expansion into (3.1). Then we omit all terms proportional to C/Oj_, where C stands for any expression with a value 
much smaller than the field strength B. We thus obtain the strong-field approximation for the first-order binding 
energy (O) 



e 



:i) 

O ,0 II 




Determining Q j_, Qu by minimization, we obtain 

2 



Q ± ^B, 



(3) 
II 



lnB - 21nhiB 



2a 
ffiB 



InB 



0(ln~ 3 B) 



(3.6) 



(3.7) 



with abbreviations a = 2 — In 2 w 1.307 and 6 = ln(7r/2) — 2 « —1.548. Thus the optimized binding energy can be 
written up to the order ln~ 2 B: 



r(l) 



(B) = i 



ln 2 B 



4 In S lnln B + 4 ln 2 ln B - Ah lnln B + 2(6 + 2) In B + b 2 



InB 



[81n 2 lnB-861nlnB + 26 2 



+C(ln- 2 B). (3.8) 

Note that the prefactor 1/n of the leading ln 2 B-term differs from a value 1/2 obtained by Landau and Lifschitz 0. 
Our value is a consequence of the harmonic trial system. The calculation of higher orders in variational perturbation 
theory should drive our value towar ds 1 /2. 

The convergence of the expansion ( |3.8| ) is quite slow. At a magnetic field strength B = 10 5 Bo, which corresponds to 
2.35 x 10 10 T = 2.35 x 10 14 G, the contribution from the first six terms is 22.87 [2 Ryd]. The next three terms suppressed 
by a factor ln _1 B contribute —2.29 [2 Ryd], while an estimate for the ln _2 B-terms yields nearly —0.3 [2 Ryd]. Thus 
we find eW(10 5 ) = 20.58 ±0.3 [2 Ryd]. This is in very good agreement with the value 20.60 [2 Ryd] obtained from the 
full treatment described in Sec. 



Ill 



Table li sts the values of the first six terms of Eq. (3.8). This shows in particular the significance of the second 
term in (3.8), which is of the same order of the leading first term, but with an opposite sign. In Fig. ||, we have 
plotted the expression £l(B) = (1/2) ln 2 B of Landau and Lifschitz 0] to illustrate that it gives far too large binding 



energies even at very large magnetic fields, e.g. at 2000Bo oc 10 12 G. Obviously, the nonleading terms in Eq. (3.8) give 
important contributions to the asymptotic behavior even at such large magnetic fields. As an peculiar property of the 
asym ptotic behavior, the absolute value of the difference between the Landau-expression £l(B) and our approximation 
(3.8) diverges with increasing magnetic field strengths B. Only the relative difference decreases. 
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IV. SUMMARY 



We have calculated the effective classical potential for the hydrogen atom in constant magnetic field, which governs 
the statistical mechanics of the system at all temperatures. At zero temperature, we find a rather accurate ground 
state energy which interpolates very well between weak and strong fields. 
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TABLE I. Perturbation coefficients up to order B for the weak-field expansions of the variational parameters and the 
binding energy in comparison to the exact ones of Ref . fl . 
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TABLE IL Example for the competing leading six terms in Eq. (E3) at B = 10 s B ps 2.35 x 10 14 G. 



(l/7r)ln 2 B 


-(4/7r)lnB lnlnB 


(4/tt) ln 2 lnB 


-(4fe/7r)lnlnB 


[2(6 + 2) /tt] InB 


b 2 /7Y 


42.1912 


-35.8181 


7.6019 


4.8173 


3.3098 


0.7632 
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